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Abstract 

An averaged version of the Lorentz force differential equation is discussed. Once the 
geometric setting has been introduced, we show that the Lorentz force equation is equivalent 
to the auto-parallel equation < L V >± x — of a projective linear connection L V on the pull- 
back vector bundle 7r*TM. Using a geometric averaging procedure, we obtain the associated 
averaged connection < L V > and we consider its auto-parallel equation < L V > ~ x — 0. We 
prove that in the ultra-relativistic limit and for narrow one-particle probability distribution 
functions whose support is invariant under the flow of the Lorentz 's force equation, the 
auto-parallel curves of the averaged connection < L V > remain close to the auto-parallel 
curves of L V. We discuss some applications of these result. The relation of the present 
work with other geometric approaches to the classical electrodynamics of point particles is 
briefly discussed. Finally, in relation with the geometric formulation of the Lorentz force, we 
introduce the notion of almost- connection, which is a non-trivial generalization of the notion 
of connection. 
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1 Introduction 

In this paper several aspects of two different topics are discussed. The first one is the notion 
of Randers space in the category of metrics with indefinite signature. The second topic is a 
geometric interpretation of the Lorentz force law and the associated averaged Lorentz dynamics. 
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Both topics are related. Indeed, one can better discuss the second one after the notion of 
semi-Randers space has been established. This is the scheme that we follow in this work. 

The notion of Randers space can be traced back to the original work of G. Randers [1]. He 
investigated an example of asymmetric space-time structure capable of incorporating the asym- 
metry of time in a geometric way. One of the consequences was a unifying theory of gravity 
and electromagnetism of charged point particles. However, the non-degeneracy of the associated 
metric tensor was not discussed in detail. This has been one problem in the mathematical defi- 
nition of a Randers space, in the case of indefinite signatures. Indeed, while for positive definite 
Randers metrics there is a satisfactory treatment (for instance [2, chapter 11]) for indefinite 
signatures, the theory of semi-Randers spaces is not universally accepted and several definitions 
are currently being used in the literature. 

There is the problem of defining a Randers-type structure which is invariant under gauge 
transformations. This is a natural requirement from a physical point of view. However, the 
combination of this with the non-degeneracy criterion of the fundamental tensor g is non-trivial. 
Combining both has lead us to define Randers spaces in the context of sheaf and pre-sheaf 
theory. The reasons to use this general framework are: 

1. To use an established mathematical formalism which is very natural to be used in dealing 
local constructions. This theory has been of fundamental importance in algebraic geometry 
[41] and in the theory of complex manifolds [42]. Grothendieck formulated a theory of fiber 
bundles with a sheaf structure [43]. 

2. One can make contact with other mathematical disciplines. In particular, one expects 
some connection between Randers-type spaces and some algebraic varieties. This point is 
work in progress. 

There are two general formalisms for indefinite Finsler spaces (that we call semi-Finsler struc- 
tures), Asanov's formalism [3] and Beem's formalism [4, 5]. In this paper, we argue why both 
frameworks and the corresponding physical interpretations are unsatisfactory, in particular when 
we try to apply them to Randers-type spaces. The main problem with Asanov's definition arises 
when one considers gauge invariance issues in Randers-type metrics; the main problem with 
Beem's formalism is that there is not known a natural definition of Randers-type metric in that 
formalism with a natural physical interpretation. 

In section 3 we provide a definition of semi-Randers space which is gauge invariant [6]. It 
has the advantage that it is taken directly from the Lorentz force equation. However, this 
definition does not correspond to a Finsler or Lagrange structure. This interpretation is natural 
in the framework for Randers-type space. All relevant geometric data is extracted from the 
semi-Riemannian metric n and from the Lorentz force equation, which in an arbitrary local 
coordinate system reads: 

d 2 a l da j da k da k I .da da. n . . . n 1 n „ ,., lN 

+ — ^vi-,-^) = 0, ,,,,* = 0,1,2,3, (1.1) 

where a : I — > M is a solution curve for r G I, V T % jk are the coefficients of the Levi-Civita 
connection W of rj and dA is the exterior derivative of the 1-form A. 

Remark. From the differential equations (1.1) one does not have enough information to 
extract a connection. Indeed, the natural object which is extracted is what we have called 
almost- connection. Nevertheless this subtlety, we use through almost all the paper the name 
connection (strictly speaking projective connection), since most of the calculations that we 
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perform are also suitable for connections (or projective connections). We formalize this notion 
later in section 11. 

In section 4 we introduce the notions of non-Linear connection associated with a second order 
differential equation. 

Section 5 is devoted to introduce a geometric formulation of the Lorentz force equation. 

Sections 6 and 7 are devoted to introduce the averaging mechanism. We have follow reference 
[7, section 4\- 

In section 8 we apply the averaging procedure to the Lorentz force equation. 

In section 9 we compare the difference of the solutions of the original system (1.1) and those 
of the auto-parallel curves of the averaged connection (starting with the same initial conditions). 
The result is that for the same initial conditions, in the ultra-relativistic limit and for narrow 
1-particle probability distribution functions, the solutions of both differential equations remain 
near each other, even after a long time evolution. Therefore the original Lorentz dynamics can 
be approximated by the averaged dynamics. 

In section 10 the limits of applicability of the averaged Lorentz model as an approximation 
of the original Lorentz force equation is discussed. Some applications of this approximation are 
also briefly mentioned. 

1. There are applications of the methods introduced here in fluid modeling of plasmas [9]. 
These applications can be of interest in accelerator physics [10] and for the description of 
the dynamics of non-neutral plasmas using fluid models [11]. 

2. The properties of a bunch of particles in an accelerator machine are measured indirectly. 
The outcomes of the measurements are associated with an averaged description of the 
bunch. Therefore, to have a theoretical model for the time evolution of these averaged 
quantities is relevant, since they are observable quantities. In some circumstances, for the 
averaged energy and momentum this is achieved by the averaged Lorentz force equation 
that we present in this work. In particular we also note that although the solutions of the 
averaged dynamics are not physical (they do not corresponds to trajectories on the unit 
hyperboloid bundle over the space-time manifold M), they can be used as a definition of 
the reference trajectory in beam dynamics [10]. 

Finally, Section 11 is devoted to formalization of the notion of almost- connection. We include 
this to complete the present work. 

Remark. In this work we are considering the problem of the mathematical description of the 
evolution of charged particles in an external Maxwell field. Therefore we are considering that 
the Maxwell equations and the Lorentz force equation hold. This differs from other geometric 
approaches to electromagnetism, where the whole theory electromagnetic theory is reconsidered 
and reformulated in a geometric way [16]. Our point of view is more conservative. Indeed, we 
have rewritten the standard theory in such a way that the averaging procedure can be applied. 
It is a surprise that in some situations, one can substitute the original Lorentz force equation 
by an averaged version. This has potential interest for several applications in beam dynamics. 
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2 Criticism of the notion of Semi-Randers Space as Space-Time 
Structure 



2.1 Randers Spaces as Space-Time Structures 

Before moving to the more specific problem of defining semi-Randers spaces, let us discuss the 
notion of semi-Finsler structure. Let M be a C°° n-dimensional manifold, TM its tangent bundle 
manifold with TM D N with projection it : N — > M, the restriction to N of the canonical 
projection tt : TM — > M. Therefore (N, it) is sub-bundle of TM. 

Let us consider the following two standard definitions of semi-Finsler structures currently 
being used in the literature: 

1. Asanov's definition [3], 

Definition 2.1 A semi-Finsler structure F defined on the n- dimensional manifold M is 
a positive, real function F : N — >]0, oo[ such that: 

(a) It is smooth in N, 

(b) It is positive homogeneous of degree 1 in y, 

F{x,Xy) = XF(x,y), VA > 0, 

(c) The vertical Hessian matrix 

l d 2 F 2 (x,y) 

is non- degenerate over N. 

9ij(x,y) is the matrix of the fundamental tensor. The set N x is called the admissible set 
of tangent vectors at x; the disjoint union N = N is the admissible set of vectors 

over M. 

In the particular case when the manifold is 4-dimensional and (gij) has signature (+, — , — , — ), 
the pair (M, F) is called Finslerian space-time. 

Remark . The strong convexity condition is fundamental for comparison results, as well 
as for the construction of connections [24, 45]. 

2. Beem's definition [4,5] 

Definition 2.2 A semi-Finsler structure defined on the n-dimensional manifold M is a 
real function L : TM — > R such that 

(a) It is smooth in the slit tangent bundle N := TM \ {0} 

(b) It is positive homogeneous of degree 2 in y, 

L(x,Xy) = X 2 L(x,y), VAX), 

(c) The Hessian matrix 



l d 2 L(x,y) 



is non- degenerate on N. 



In the particular case when the manifold is 4-dimensional and gij has signature (+,—,—,—), the 
pair (M, L) is a Finslerian space-time. 
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2.2 Comparison of Asanov's and Beems's Formalism 

Some differences between the above definitions are highlighted below: 

1. In Beem's framework there is a geometric definition of light-like vectors and it is possible 
to derive Finslerian geodesies, including light-like geodesies, from a variational principle 
[12]. By construction, in Asanov's formalism it is not possible to do that in an invariant 
way, because light-like vectors are excluded in the formalism from the beginning, since 
nothing is said about how to extend the function F 2 from N to TM. 

2. In Asanov's framework, given a parameterized curve a : I — > M, I C Ron the semi- 
Finsler manifold (M, F) such that & £ N for all t El, the length functional acting on a is 
given by the line integral 

f'tmax 

a(t) — > £ A (a) := F(a(t),a(t))dt, I =}t min , t max [. (2.3) 

Due to the homogeneity condition of the Finsler function F, £a(o~) is a re-parametrization 
invariant functional. On the other hand, if we consider Beem's definition, the energy 
functional is given by the following expression [12]: 

a(t) — > £ B (a) := L(a(t),a(t))dt, I =}t min , t max [. (2.4) 

Formulated in this way, Beem's energy functional is not re-parametrization invariant, 
because the fundamental function L is homogeneous of degree two in y. 

3. A third difference emerges when we consider the category of Randers-type spaces: 
Definition 2.3 (semi-Randers Space as semi- Finsler Space) 

In Asanov's framework, a semi-Randers space is characterized by a semi- Finsler function 
of the form: 

F(x, y) = \jr){x)ijtyi + A(x, y), (2.5) 

where rjij(x)dx t <8> dx 3 is a semi-Riemannian metric defined on M and A(x,y) := Ai(x)y l 
is the value of the action of the 1-form A(x) = Aidx 1 acting on y £ T^M. 

In the positive definite case and when rj is a Riemannian metric, the requirement that (gij) 
is not non-degenerate implies that the 1-form (A\, ...,A n ) is bounded by rj: 

AiAjrf 3 < 1, rf k r] kj = &j. 

The indefinite case is quite different, since there is not a natural Riemannian metric that in- 
duces a norm in the space of homomorphisms. Therefore, the criterion for non-degeneracy 
becomes unclear. 

Secondly, for both the positive and indefinite metric, only the variation of the length func- 
tional (2.3) (not directly the integrand itself) is invariant under the gauge transformation 
A — > A + dX; the Finsler function (2.5) is not gauge invariant as well. Even in the case 
that we could define a metric and norm, transforming the 1-form A by a gauge transfor- 
mation can change the norm and therefore the hessian (g^) can become degenerated. This 
fact suggests that while the variation of the length functional has a physical meaning, the 
fundamental function F does not allow this interpretation. 
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On the other hand, the notion of Randers space in Beem's formalism is even more prob- 
lematic. In this case there is not a formulation of semi-Randers space (because eq. (2.5) is 
of first order in y). This suggests that a proper formulation of the notion of semi-Randers 
space requires going beyond metric structures. 

4. It is interesting to have a definition of semi-Randers structure capable of taking into 
account light-like trajectories for charged particles. As we have seen, Asanov's treatment 
is not able to consider light-like vectors. On the other hand, Beem's formalism is not 
capable of considering this problem, since there is not a known Randers-type structure 
in Beem's formalism. However, the asymptotic expansion of the ultra-relativistic charged 
cold fluid model presented in [13] is an example where this light-like trajectories appear 
naturally. In that model, the leading order contribution to the main velocity field of the 
charged cold fluid model is a light-like velocity vector field, interacting with the external 
electromagnetic field; perturbative corrections change the velocity vector field to a time-like 
vector field. 

The above observations make it reasonable to introduce a non-metric interpretation for semi- 
Randers spaces. The option that we have adopted has been to formalize a geometric struc- 
ture from the geometric and physical data that we have: the Lorentz force equation and the 
Lorentzian metric n. 

3 Non-Lagrangian Notion of Semi-Randers Space 
3.1 Non-Lagrangian Notion of Semi-Randers Space 

Given a manifold M, let us assume the existence of a smooth semi-Riemannian structure rj on 
M. This implies that the function 

■n : TM x TM — > R 

(X, Y) ' V THjix^Yi, X, Y e T X M 

is smooth on the variables x, X l ,Y J . Since we will use the square root \/rjij(x) X l Yi , we also 
require that: 

1. y/rnj(x) XiYj is smooth in "N := LUtm {v G T x M, mj (x) X l Y> > 0}. 

2. v /-r /lJ {x)X i Y^ is smooth in "N := \J xeTM iv G T X M, mj (x) X^ < }. 

The null-cone is ^NC := U^eM iv ^ T^M | n(y, y) =0}. Let us assume the existence of a closed 
2-form F on M. We propose a notion of semi-Randers space based on the following 

Definition 3.1 A semi-Randers space consists of a triplet (M, rj, F), where M is a space-time 
manifold, rj is a semi-Riemannian metric continuous on M and smooth on TM \ ^NC := 
r 'N U and a 2-form F G f\ 2 M such that dF = 0. 

F is on the second de Rham cohomology group H 2 (M.). Due to Poincare's lemma, there is a 
locally smooth 1-form A such that dA = F. Any pair of locally smooth 1-forms A and A such 
that A = A + dX, with A a locally smooth real function defined on the given open neighborhood, 
are said to be equivalent and produce under exterior derivative the same cohomology class 
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[F] G H 2 (M.) that contains the element F: d(dX + A) = dA = F. Note that we are speaking of 
locally smooth 1-forms A and of globally smooth 2- forms F. 

Therefore, instead of giving F, one can consider the equivalence class of 1-forms: 

[A] := {A = A + dX, dA = F in the intersection of the opens sets where A and A are defined}, 

with A a locally smooth 1-form defined on an open set U C M. Then, if two 1-forms ^ A and 
v A are in [A], d^A = d v A, where [i and v are in a family of an arbitrary family of indices. 
This implies that d( ^ A — v A) = 0. For each point x of the open neighborhood where both 
potentials are defined, there is a locally smooth function defined in an open neighborhood of x 
XJ(x) C such that (^A - v A) = d(^ v X) (using Poincare's lemma again). 

Based on these arguments we give an alternative definition of semi-Randers space: 

Definition 3.2 A semi-Randers space consists of a triplet (M, rj, [A]), where M is a space-time 
manifold, rj is a semi-Riemannian metric continuous on M and smooth on TM \ LC and the 
class of locally smooth 1-forms A G [A] is defined such that dA = F for any A G F. 

Proposition 3.3 Both definitions of semi-Randers space are equivalent. 

Proof. We showed already one of the directions of the equivalence. To show the other direction, 
one needs to construct locally 1-forms which produce the required 2-form F. This is achieve by 
the Poincare lemma in a start shape domain [27, pg 155-156]: 

fl n-l 

A(x) = ( tJ2 xkF kj(tx)dt)dx j . 

□ 

We adopt definition 3.2, since it has the advantage that it allows us to discuss some local 
issues related with the inverse variational problem of the Lorentz force equation. Essentially, 
this is the reason that even if the topology of M is trivial, the 1-forms A are only locally smooth. 
For two neighborhoods M U and ^U, the potentials are ^A and v A. Since they are in the same 
equivalence class [A] in the intersection U XJ n / they are related by ^A = d(^ u X) + U A, 
where >* U X is also a locally smooth function. 

Remark. Note that this construction does not imply any additional restriction on the group 
fl^M). 



Let us denote the pre-sheaf of locally smooth functions over M by /\f oc M. For each of 
the representatives A £ [A] G /\j oc M. there is on M a function Fa defined by the following 
expression: 

F , \ = / y / Vij(x)y i y j + Aj(x)y* for r)ij{x)y l yi > 0, 
M ,!J> I fW + AWl/' for %(x)yV<0. 1 • J 

The following properties follow easily from definition (3.1) and from the definition Fa'- 

Proposition 3.4 Let (M,t/, [A]) be a semi-Randers space with rj a semi-Riemannian metric, 
A G [A] and Fa given by equation (3.1). Then, 

1. On the null cone NC^ := {y G T X M | % (x) y l yi = 0} Fa is of class C° , for rj G C° and 
A G AL M - 
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2. The subset where r) x (y, y) / is an open subset of T X M. and Fa is smooth on T X M\NC X , 
for r] smooth and A G Aj oc M. 

3. There is no constraint on the non- degeneracy of the fundamental tensor gij. Therefore, it 
is not required that any representative A of [A] be bounded by 1. 

4- The function Fa is positive homogeneous of degree 1 in y. 

Remark. Without the requirement that the fundamental tensor g^ be non-degenerate, it is not 
possible to define a geodesic equation from a variational principle [14]. 

3.2 Variational Principle Based on semi-Randers Spaces 

Let Q(M) be the set of smooth curves on M with time-like tangent vector field. The functional 
acting on a is given by the integral: 

£ Fa — > R 

f F A (a(r),a{T))dT (3.2) 

with r the proper time associated with r\ along the curve a. The functional is gauge invariant 
up to a constant: if we choose another gauge potential A = A + dX, then £f a {o~) = £f a (&) + 
constant, the constant coming from the boundary terms of the integral. Therefore, the variation 
of the functional is well defined on a given semi-Randers space (M, i], [A]). We denote the 
corresponding functional acting on a by £f(o~). 

Remark. The constant can be arbitrary large. Also note that the entire variation of the 
curve of the geodesic a must be inside the domain of a given constant. 

The condition that guarantees the construction of the first variation formula and the existence 
and uniqueness of the corresponding geodesies is that the vertical Hessian gij must be non- 
degenerate [14]. However, given a particular potential A £ [A] one can not guarantee that the 
Hessian of Fa is non-degenerate. However, due to the possibility of doing gauge transformations 
in the gauge potential A(x) i-> A(x) + d\(x) we have: 

Proposition 3.5 Let (M, n, [A]) be a semi-Randers space. Assume that the image of each 
solution a on the manifold M is a compact subset. Then 

1. There is a representative A 6 [A] such that the Hessian of the functional Fa is non- 
degenerate. 

2. The functional (3.2) is well defined on the Randers space (M,r/, [A]) ; except for a constant 
depending on the representative A £ [A]. 

3. If the geodesic curves are parameterized by the proper time associated with the Lorentzian 
metric 7], the Euler- Lagrange equations of the functional Fa are the Lorentz force equations. 

Proof: 

1. Using the gauge invariance of £f(o~) up to a constant, we can obtain locally an element 
A G [A] such that AiAj rf 3 < 1 in an open neighborhood in the following way. Consider 
that we start with a 1-form A which is not bounded by 1. The 1-form A[x) = A{x) + d\{x) 
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is also a representative of [A]. The requitement that the hessian of A is non-degenerate is, 
using a generalization of the condition [2, pg 289]: 



< 



2 ^>)y l + vWF Ai(s)Aj(a;)) 



On the unit tangent hyperboloid S z this condition reads: 



< 



2 + 



i+ Ai(x)y % 

Let us assume that (if it is negative, the treatment is similar) that 

In order to be able to write down this condition one needs that 1 + Ai(x)y % / 0; the region 
where this does not hold is the intersection of the hyperplane 

P., :={yG T X W | 1+ A^y 1 = } n H x 

with the unit hyperboloid S x . The intersection is such that P X C\ T* x C {y G T* x \ F(x, y) = 
0}. 

Let us write in detail the above condition of positiveness (Ai = Ai(x) + diX(x)): 
^ _ n , (Ai{x) + di\(x))y l + rf' (A^x) + ^A(x)) (^(x) + djX(x)) 



e z (x,y) = 2 + 



1+ (^(ar) +d i \{x))y i 



If y G S x , then for /? 7^ 1, fly is not on the unit hyperboloid. The situation is different 
if y G NC X . Then, /3y G NC X . Now we make the approximation of ~E X — > NC^ in 
the asymptotic limit y° — > 00. In this approximation, one can perform limits in the 
expression for e 2 (x,y). In particular, one can consider e 2 (x,f3y) for large f3: 

lim e 2 (x By) = 2+ lim ^ + d jMtK + ^ M g ) + ^)) Wf) + ^ A(x) ) - 
/?->oo 1 /3^oo 1+ (Ai(x) +d i \(x))py i 

= 2+ (A t (x)+d t \(x))W 3 

^00 1+ (^i(x) + <%A(x))/V 

That this limit is well defined for large enough y° can be seen in the following way. Let us 
assume that 1 + (Ai(x) + di\{x))y l = 0. Then, we consider the expression 1 + (Ai(x) + 
di\(x))(3y l for j3 » 1. It is impossible that the second expression is zero except if (3 = 1 
and 1 + (Ai(x) + d i X(x))py i = 

One should prove that e 2 (x,y) is positive for any y G ~S X . This is achieve because e 2 (x,y) 
is invariant. Therefore, we can change to a coordinate system where y° is arbitrary and 
the value of the bound does not change. 

2. When the curve a whose image is defined over several open sets, for instance and 
it is evaluated using several local potentials. If they are representatives p 'A u A and there 
is a contribution coming from boundary terms, due to the evaluation of the d\ terms on 
some boundary points in the intersection ^Un This differences does not contribute to 
the first variation of the functional. Therefore, the first variation of (3.2) exists and does 
not depends on the representative. The corresponding extremal curves exist and they are 
unique. They correspond to the solutions of the Lorentz force equation. 
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3. If the image is a compact set, the number of sets M U that we need is finite. Then we 
obtain a globally defined section of Aj M. 

4. Once a global locally differentiable 1-form with the required properties is obtained globally 
over the variation Var(a) of a : I — > M, one can follow the standard proof of the 
deduction of the Lorentz force equation from the variation of a functional [1], [36, pg 47- 
52]. However, one of the formal differences with standard approaches is that we restrict 
the variational approach to compact regions. □ 

Remarks 

1. It is important to notice that the non-degeneracy of the metric gij not necessarily implies 
that it has the same signature than the semi-Riemannian metric rj. Further investigations 
are required to determine the criteria for conservation. 

2. If the curve a is parameterized with respect to a parameter such that the Finslerian arc- 
length F(a, &) is constant along the geodesic, the geodesic equations have a complicated 
form (see for instance [2, pg. 296]) and they are not invariant under arbitrary gauge 
transformations of the 1-form A — > A + dX, X € Fi oc {M). 

3. The fact that we are speaking of local data makes natural to consider notions from pre- 
sheaf theory [27 chapter 6] as the basic ingredient in the definition of Randers spaces. 
Sheaf theory is a theory that allows to treat local objects defined on sheafs, for example 
germs of locally smooth functions or locally smooth differential forms [28, chapter 2}. In 
this case, the definition of the class [A] is on the second class of the sheaf cohomology on 
the manifold M. 

With definition (3.2) at hand the problem of how to introduce the gauge symmetry in a 
Randers geometry is solved. The price to pay is: 

1. The class [A] and the Riemannian metric rj are unrelated geometric objects. This is in 
contradiction with the spirit of Randers spaces as a space-time asymmetric structure. 

2. One needs to consider locally smooth sections instead of globally defined sections. This 
happens even if the topology is trivial. 

We have seen that given a 1-form A 6 [A] one has to work hard to find another representative 
A £ [A] such that Douglas's theorem [14] holds almost 

On the other hand, light-like trajectories can not be considered in this formalism for semi- 
Randers spaces in a natural way. Beem's formalism allows the treatment of light-like geodesies 
as extremal curves of an energy functional for some examples of indefinite Finsler space-times. 
However, it is not known a semi-Randers functional in Beem's formalism. 

4 Non-Linear Connection Associated with a Second Order Dif- 
ferential Equation 

Let M be an n-dimensional smooth manifold. A natural coordinate system on the tangent 
bundle ir : TM — > M is constructed in the following way. Let (x, U) be a local coordinate 
system on M , where U C M is an open sub-set M and x : U — > R n a local coordinate system. 
An arbitrary tangent vector at the point p £ U is of the form X p := X = X k -^\ p . The 
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local coordinates associated with the tangent vector X p 6 T Z M C TM are (x k ,y k ). From the 
imbedding e : N i— > TM e(N) acquires the induced differential structure from TM; e(N) is 
denote by N. Indeed, N is a sub-bundle of the tangent bundle TM of maximal rank. 

We recall the following notion of connection [15, pg 314]. Let tt : N — > M the differential 
function dir : TN — > TM. Then the vertical bundle is V = kerdn. 

Definition 4.1 A connection in the sense of Ehresmann is a distribution H C TN such that: 

1. There is a decomposition at each point u G N, T U N = jH u © V u . 

2. The horizontal lift exists for any curve t i-> a{t) £ M, t\ < t < t2 and is defined for each 
t G R. 



Let us consider a set of n = dim{NL) second order differential equations. The solutions are 
parameterized curves on M. Assume that the system of differential equations describes the flow 
of a vector field G x S TTN. In particular, the system of differential equations has the following 
form: 

^--G i (x,^) = 0, i,j = l,...,n, (4.1) 

for a given parameter t. This system of differential equations is equivalent to a first order system 
of differential equations on N, 

{ Hr = y 3 ^ i,J = l,...,n. 

The coefficients G l (x, y) are called spray coefficients if they are homogeneous functions of degree 
one on the coordinate y; the general case where they are not homogeneous are called semi-spray 
coefficients. G 1 (x,y) are such that they transform under a change of natural local coordinates 
on N, induced from changes of coordinates on M, in such a way that the system of differen- 
tial equations (4.1) remains invariant in form. Explicitly, if the change in the local natural 
coordinates on the manifold N is of the form: 

x 1 = x l (x), 

f = By'- 

The associated co-frame transforms in the following way: 

/ <'■<" = 

The induced transformation in the associated system of differential equations are 

?£-G\x,y) = ?£-G i (x,y) = 0, 
which implies the transformation of the coefficients G l (x, y), 

j,k 3 
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The vertical distribution V admits a local holonomic basis given by 



{^T,-..,^}, '••/ 1 "• (4-3) 

Using these spray coefficients it is possible to define a horizontal n-dimensional distribution of 
the fiber bundle TN — > N, 

{— — } -i_-=A_^lA i - = i n (4 4) 

such that it generates a supplementary distribution to the vertical distribution. The non-linear 
connection coefficients N l k{x,y) are defined by the relation 

G l (x,y) :=y k N l k (x,y). 

For a spray, the connection coefficients of the non-linear connection are: 



N\(x,y) 



dy k 



Since the spray coefficients G % transform in a well defined way, the non-linear connection coeffi- 
cients N l j(x,y) also transform in a characteristic form [2, 21]: 

(■<" = v l = %v 3 ) ■ K>(*>v) = N T^y) f§ (x) - 

Given a spray G l (x,y), the non-zero coefficients correspond to the covariant derivative of 
horizontal sections of TTN along horizontal sections of TN and they are given by the Hessian 
of the spray: 

x F i. k ._ 1 d 2 G\x,y) _ 
jfe 2 dyidy k 

All the other coefficients are zero. This type of connection resembles the Berwald connection 
used in Finsler geometry. From the geometry of sprays, it is a natural connection. Note that, 
while x r* fc can be associated with a linear connection in a given vector bundle, the connection 
coefficients N % j(x,y) are not (this is why they are called non-linear connections). 

Given the non-linear connection, one can define the horizontal lift of the tangent vectors; the 
horizontal lift of X = X l di G T^M to the space T M N is defined to be h{X) = X 1 -^. This 
lifting is defined here using local coordinates. However, an intrinsic definition can be found 
in [40]. From this basic definition, one can define the horizontal lift of vector sections and in 
general, tensor sections of the corresponding bundles. 



The pull-back 7r*TM — > N of the bundle TM is the minimal sub-manifold of the cartesian 
product N x TM such that the following equivalence relation holds: for every u G N and 
(u,£) G 7r ] ^ 1 (u), (u,£) G 7r*TM iff 7r o tt2(u, £) = 7r(u); the projection on the first and second 
factors are denoted by 

TTi : 7r*TM — ► N, (u, — > u, 
vr 2 : 7r*TM — »• TM, (u, — >• 
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The pull-back bundle 7r*TM — > N is such that the following diagram commutes, 



7T*TM ■ 



N- 



7r*TM — > N is a real vector bundle, with fibers homeomorphic to T X M.. Similar pull-back 
bundles can be constructed from other tensor bundles over M, for instance 7r*T*M — > N. 

Given the non- linear connection on the bundle TN — > N, there are several related linear 
connections on the pull-back bundle 7r*TM — > N. All these connection coefficients are based on 
the coefficients* r* j k . In particular we can stipulate the following connection on 7r*TM: 

V^vr*Z := x T(x, y y jk Z j ir* ei , V v ir*Z := 0, V G V; (4.5) 

here {n*ei, i = 0, n — 1 } is a local frame for sections r(7r*TM) The can be generalized this 
construction to generic tensor bundles over M. 



5 Geometric Formulation of the Lorentz Force Equation 

5.1 The non-linear Connection Associated with the Lorentz Force Equation 

Let us consider a semi-Randers space (M,?7, [A]), and the bundle N — > M, N := U^gM {y G 
T X ~M, rj(y, y) > 0} C TM. Then, the following diagram commutes: 

TM 



N— ^M. 

where i is the following natural embedding 

l : £ — > TM 

(x,y) ^ (x,y), x G M, y € ~E X . 

ir is the restriction of 7r to the sub-bundle N. Let us consider the differential map dir : TN — > 
TM. Recall that the vertical bundle is defined as the kernel V := ker(d7r); at each point one 
has ker(d,Tt\ u ) := V u , with u G N. 

The system of second order differential equations (1.1) determines a special type of vector 
field on N called spray. It is well known that a spray defines an Ehresmann connection [16]. Let 
us denote by 7](Z,Y) := rjij(x)Z l Y J . 

Definition 5.1 Let (M.,i],[A]) be a semi-Randers space. For each tangent vector y £ T X M 
with r](y, y) > 0, there are defined the following functions: 

L V jk (x,y) = T% fc + - * (F l J (x)y m Vmk + F i k (x)y m 7 lmj )+ 

2Vv{y,y) 

y m 1 

+F '"' (I) v^ ( ''''"^' , '* ! '* ! ' ) ' (5 - 1) 

where r/(y,y) is a short way to write r]ij(x)y t y 3 , j k , (i,j,k = 0, 1,2, ...,n) are the connection 
coefficients of the Levi-Civita connection T 'V in a local frame, F^- := diAj — djA-i and F l j = 
rf k F k j, for any representative A G [A]. 
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The unit hyperboloid sub-bundle XI acquires an induced connection, which connection coeffi- 
cients are 

L r ifc (x,y)| s = ^T i jk + ^(F i ] (x)y m 7 ]mk + F i k (x)y m 7 ]mj )+ 

y m 

+ F 'mWyfe - VjsVkW s y 1 )- 

The structure of this functions is clear: T* fc are the connection coefficients of the Lorentzian 
metric 77; the other two terms are tensorial. Indeed, one can define the following expressions: 

V jk = - 2- _ (F i jy m Vmk + F l k y m r, mj ), 

v m ( 1 \ 

Therefore, 

L-pi VT^i 1 r f^ 1 

1 jk — 1 jfe ~r J- jk ~t~ J-i jk: 

With the functions U- k and Tj fe , one can construct the following maps: 

L u : T U N x T M N — > T U N 

(X,Y)^ L i jk [x,y)XiY h -^- i . 
The second operator that we define is: 

T u : T U N x T M N — > T M N 

■u = (x, y) and X, Y are arbitrary tangent vectors X, Y G T U N. This can be generalized to 
homomorphisms acting on vector sections. 

Note the following elementary property: 

T U (Y, Y) = 0, Vy G T X M, Y = y^, « = (x, y). 
However, T u (-, Y) / and that T U (Z, Z) / for arbitrary Z. 



5.2 The Koszul Connection on TN Associated with the Lorentz Force Equa- 
tion 

Let {ei, e n } be a local basis for the sections of the frame bundle associated with the tangent 
bundle TTM and let us assume that each a is a time-like tangent vector at x (therefore, the met- 
ric cannot be diagonal in this basis, since n is a semi-Riemannian metric). Then {ir*ei, ...,ir*e n } 
is a local frame for the fiber ir^ 1 C 7r*TM, u G N; {hi, h n } is the local frame of the horizontal 
distribution % u C T U N obtained by horizontal lift h\ = h(ei) and {v±, v n } is a local frame 
for the vertical distribution V u C T U N. 

Given the set of functions { L T l j k , i,j, k = 0, ■ n— 1} and the non-linear connection associated 
with the system of differential equations, there is an associated Koszul connection on TN [17]. 

Proposition 5.2 Let (M, 77, [A]) be a semi-Randers space and M be of dimension n. There is 
defined a covariant derivative L D on TN determined by the following conditions: 



15 



1. For each X G % u and Z G VH U 

L D X Z = X kL V jk {x,y)Zi h u X = X 1 hi\ u , Z = Z t h i \ v , 

with {hi} a local frame for the horizontal distribution, u = (x,y) and v an arbitrary point 
of an open set OcM containing u. 

2. The covariant derivative of arbitrary sections Z G TTN along vertical direction is zero: 

L D v Z = 0, VVgV, ZgTTN. 

3. The covariant derivative L D is symmetric, i.e, has zero horizontal torsion: 

L D V V- L D V U- [U,V] = 0, VU,V eH. 

4. For all X eH u and Z G TV U , L D X Z = 0. 

Proof: Let us consider the Finsler geodesic equation associated with the semi-Randers space 
Fa = \/i]ij(x)y t y : > + Ai(x)yi but parameterized using the arc-length of the Lorentzian metric n: 



d 2 x % „ dxi dx k a, ias I A x dx.dx k 

where F = dA is the exterior differential of the 1-form A and r](X, Z) = r]ij(x)X l ZK From these 
equations, we can read the value of the semi-spray coefficients [34]: 



G\x,y) = jk (x)yiy k + V ^ (dA) jk (x)^y~^y k . 
Taking the first and second derivatives with respect to y, we obtain: 

l^-G l (x, y) = y l + rf\x) (dA) lm (x) ^L^ Vjs (x)y s y m + rf\x){dA) lj {x)^(^y). 

\w k W Gi{x ' v) = Vrk3{x) - v il (x)(dA) lm (x) 2(^)3/2 y s y m m P (x) y*>+ 

+rf\x){dA)im{x) n( r L^ r ljk (x)y™ + rf\x){dA) lk {x) l^ m s y s + 
+rf l (x) (dA)ij (a) -= =- Vks(x)y s . 

2(V'?(y>y)) 

One can check that l^-^G^x^y) = L T l jk (x, y). 

From the definition of these coefficients one can check (following a standard procedure, for 
instance in [16]) that there is a splitting of the tangent vector spaces T U N for each u G N. In 
addition, one can check that the connection coefficients are symmetric L Y i j k = L T i k j. The fact 
that the covariant derivative along the vertical directions is zero is an additional hypothesis used 
to make the covariant derivative unique. □ 

Proposition 5.3 The following properties hold: 

1. The Lorentz connection L D is invariant under gauge transformations A — > A + dX of the 
1-form A{x) = Ai(x)dx l . 
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2. Given a point x £ M, D admits a normal coordinate system centered at x which coincides 
with the normal coordinate system associated with W centered at x iffF(x) = 0. 

Proof: 

1. The first property comes from the fact that all the geometric objects appearing in the 
connection coefficients L Y % jk are gauge invariant. 

2. If L D is affine, for any point x £ M there is a coordinate system where the connection 
coefficients are zero L P j k = 0. This implies that L T l jky^y k = at the point x. Due to 
the decomposition L V* fc = v Vj k + T l j k + L l j k this is equivalent to: 

= rr i jk + T i jk + L i jk )yiy k , Vj,eS,. 

Since the transversality condition holds (T l jky^y k = 0), one obtains: 

= (' l r i jk +L i jk )yiy k , Vj/eS,. 

Assume that there is a normal coordinate system centered at x for L V and that this 
coordinate system coincides with the normal coordinate system associated with W. Then 
at the point x, one has the relation 

( L r i jk yiy k )= F'jyi = 0,VyeT x M. 

This last condition is strong enough to imply F = 0. 

□ 

Remark. In the positive case, it is well known that the requirement that the Chern connection 
of a Randers space lives on the manifold M is that the 1-form A must be parallel V VA = 0. This 
is a stronger condition than dA = 0. The parallel condition indicates that the space is Berwald; 
the closeness condition indicates that the space is Douglas [2, pg 304]; Douglas's spaces are such 
that they have the same geodesies than the underlying Riemannian metric rj. One can easily 
move the proofs to the Lorentzian and indefinite setting, if one adopt Asanov's framework. 

Corollary 5.4 Let (M, rj, [A]) be a semi-Randers space. Then the Lorentz force equation can 
be written as 

L D~x = 0, 

where x : I — > M is a time-like curve parameterized with respect to the proper time of the 
Lorentzian metric rj, x is the horizontal lift on N and L V is the non-linear connection determined 
by the system of differential equations (1.1). 

Proof: A solution of the auto-parallel condition of the Lorentz connection defines a curve on 
N given by (x,y)(t) = (x(t),x(t)). Projecting this curve into M using n, one obtains a curve 
x(t) which is a solution of the Lorentz law. □ 

5.3 Lorentz Connection on the Pull-back bundle 7r*TM 

We introduce the third framework, which will be directly used later to define the averaged con- 
nection. Given the non-linear connection L D on TN — > N, there is a natural linear connection 
on the pull-back bundle 7r*TM — > N that we denote by L V characterized by the following: 
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Proposition 5.5 The linear connection L V on the pull-back bundle 7r*TM — > N is such that: 

1. L V on 7r*TM — > N is a symmetric connection : 

L V^vr*y- L Vy7t*X -tt*[X,Y] = 0, (5.2) 

where X,Y G TTM, X,Y G TTN are horizontal lifts of X,Y G TTM to TTN, 
r/(X, X) > and r/(Y, F) > 0. 

2. The covariant derivative along vertical directions of sections of 7r*TM are zero: 

L V Vi 7r*e k =0, (j,k = l,...,n). (5.3) 



L V hi n*e k = L r) k {x,y)7r*ei, k = 1, n). (5.4) 



3. The covariant derivative along horizontal directions is given by the formula: 

■ jk 

4- By definition the covariant derivative of a function f G .F(N) is given by 

L Vxf:=X(f), VXGT U N. 

Proof: One checks by direct computation that these relations define a covariant derivative on 
7r*TM and that they are self-consistent. A general covariant derivative can be expressed in 
terms of the connection 1-forms: 

J 3 := L r Jk dx k + L T' Jk 5y\ 

where we have used a local frame of 1-forms {dx 1 , dx n , Sy 1 , Sy n }. Since the covariant 
derivative of sections on 7r*TM along vertical directions is zero, one obtains 

L T jk 5y k = L T 3k = 

at each point (x, y) G N. Since the torsion tensor is zero: 

L-pi L-pi 

jk 1 kj ■ 

Therefore, we only have to provide the rule of how to derive sections along horizontal directions. 
Since the coefficients given by formula (5.1) are symmetric, this rule is consistent with the 
torsion-free condition. Finally, we want that the covariant derivative to be a local operator. 
This is satisfied by (5.5), which guaranties that L X7 satisfies the Leibnitz rule. □ 

Corollary 5.6 Let M, N, 7r*TM and L X7 be as before. Then the auto-parallel curves of the 
linear Lorentz connection L V are in one to one correspondence with the solutions of the Lorentz 
force equation: 

L V-ir*x = L D ± ± = 0, x = ^-. 

dr 



L v 



Proof: If in some coordinate system the Lorentz connection L V has the connection coefficients 
jk'- 

(^{Lrj ,* dx ( T ) \Y f f^jr) L i , dx{T) ^ dxi{T)dx k {T) \ 

_/dV(r) / _ u 1 dx^r)dx-(r) dx>(r) 

~\ dr* V {dA)lm 2{v{ ^^ )) 3/2 7] > s dr dr Vkp dr + 
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U/ias 1 dx m (r) 1 dx s (r) 

+rf l {dA) lm Vjk —^ + rf\dA) lk ■ r?js _^+ 

Using rj( dx ^ , dx ^ ) = 1), the above expression simplifies: 



(V( 



V d£(r) 



d 2 x*(r) /__,_■ ii, 1A s 1 ^ s (r) , 7/ ,^ 1 dx s (T)\dx j (T)dx k (T)\ 



□ 



6 The Average Operator Associated with a Family of Automor- 
phisms 

6.1 Average of a Family of Automorphisms 

The averaged connection was introduced in the context of positive definite Finsler geometry in 
[7]. However, in this work we need to formulate the theory for arbitrary linear connections on 
the bundle 7r*rM -)• N, where N — ► M is a sub-bundle of the tangent bundle TM — > M. 

Let 7r*,7ri,7T2 be the canonical projections of the pull-back bundle tt*T^ p '^M -)• N, T&'^M 
being the tensor bundle of type (p,q) over M, ir^T^M the fiber over u G N, Ti P ' 9) M the 
tensor space over x G M and S x a generic element of T^'^M; S u is the evaluation of the section 
S G T^T^M) at the point u G N. 

The averaging operation requires the following structure: 

1. A family of non-intersecting, oriented sub-manifolds 

:= C T^M, x G U C M} C T^M; 

the particular case £ = E^ 1 ' ), = is relevant case for us. 

2. A measure at each point x, which is an element f{x,y)u x {y) G /\ m ^) x , where m is the 
dimension of T* x and / s := /(x, y) has compact support; supp(f x ) C X^; / x : Xl^ — > [0, oo] 
is required to have compact support. 

j](p.g) — >. M is a bundle over M; the fiber ~E X can be non-compact. However, using local bump 
functions one can reduce the domain of / to be compact. 

We need to introduce some additional notation. For each tensor S z G T^'^M and v G 7r _1 (z), 
z G U C M the following isomorphisms are defined: 

TT 2 \v ■ 7T*T^M — »• T<f'^M, S v — »• S 2 

< : T^M — > <T^)M, 5* — ► n* v S z . 
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Definition 6.1 Consider the family of automorphims, 

{A w : tt*TM — ► tt*TM, w G 7r _1 (a;)}. 
The averaged operator of this family is the automorphism 

< A > x : T X M — »• T,M 

Sx — Tr^ \ i I n 2\uA u K) ■ S x , 
vol(E x ) \ J^ x J 

u G 7T _1 (x), S x G r^M. 
The integral operation means the following: 

S^) f(x,u) dvol(x,u), 

where dvol(x,u) is the measure on the fiber 7r _1 (a;) C N and must have compact support. We 
will denote the averaged automorphisms by symbols between brackets. 



6.2 Examples of Geometric Structures which Provide an Averaging Proce- 
dure 

1. Lorentz Structure [34]. The geometric data is a Lorentzian metric rj on M. The disjoint 
union of the family of sub-manifolds S x C T X M defines the fibre bundle ir : £ — > M, 
which we call the unit hyperboloid bundle over x G M: 

£:= y{S,cT,M}, V x :={yeT x M\r)(y,y) = l} 

The manifold ~E X is non-compact and oriented. The measure on ~E X is given by the following 
(n — l)-form 

co x (y) := f(x,y)^j±dy 1 A • • • Ady n ~\ y° = y\x\x\ ...,x n ~\y\ ...,y n ~ 1 ), 

where f x (y) = f(x,y) is positive and with compact support on S^; the function y° defines 
the parameterized hypersurface S x C T^M. 

Note that the Lorentzian metric rj does not determine the manifolds {£;c, % G M}. For 
instance, one can consider £ to be the collection of null cones over M: 

NC := [J {NC, C T X M}, NC, := {y G T X M | r,(y, y) = 0}. 

7r : NC — t- M is the light-cone bundle over M and NC X is the light-cone over x; on the 
other hand, e : NC ^ TM is a sub-bundle of TM — > M. 

2. Finsler Structure [2, 32]. In this case, the Finsler function F(x,y) defines a fundamental 
tensor gij(x,y) = \^~§^r§^- which is positive definite, homogeneous of degree zero on y, 
smooth and living on the sub-bundle N := TM \ {0}. In this context, £ is defined as the 
disjoint union, 

£ := □ {I, C T-jM}, I x := {y G T X M | F(y, y) = l}. 
xeM 
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S is the indicatrix bundle over M; the manifold I x is compact and strictly convex for each 
x 6 M and is indicatrix at x. The volume form is 

dvol(x,y) := f(x, y)Jg ^ dy 1 A • • • A dy n ~\ y° = y (y\ ...,y n - 1 ). 

3. Symplectic Structures [25]. In this case, there is defined on T*M a non-degenerate, closed 
2-form oj. Due to Darboux's theorem [25, pg 246], there is a canonical local coordinate 
system of T*M such that the symplectic form oj can be written as 

n 

oj = ^2 dp 1 A dq l . 

i=l 

Associated with oj there is defined on the dual tangent bundle T*M a volume 2n-form: 

S = oj A •••Aw. 
In the canonical coordinates (q,p), the 2n-differential form is 

S(p, q) = dp 1 A • • -dp n A dq 1 A • • -dq n . 

Therefore, let us explore the following construction. Let us assume the existence of a 
non-where zero vector field V (therefore the Euler characteristic of TM must be different 
than zero). Then we can construct the following form (2n — l)-form, 

u q [p) = LyS, V G TT(T*M) 

which is a non-degenerate (2n — l)-form, whose value on V is zero. Let us choose a dis- 
tribution of commuting vector fields, {Xi, [Xi,Xj] = 0, i,j = l,...,2n — 1} locally 
complementary to V and such that ...,X2 n -i, V} is a local frame of TM. The dis- 
tribution {X±, X2n-i} is integrable and dvol(x,y) is a volume form on the integral 
manifold S^. On the other hand, it is a fibered manifold: 

tt:S ± := □ S\ — > M. 
xeM 

The interesting thing about this example is that we can only construct local averaging 
procedures; the overlapping of the averaging operation in different open sets is non-trivial 
and in general one needs more structure to define consistently the averaging procedure 
globally. 

4. Hermitian Vector Bundles [35]. The construction is similar to the one in the Finslerian 
case. The sub-manifolds ~S X are defined as: 

V x = {yeT x M\H{y,y) = l}, 

where H is the hermitian structure on M. Therefore 

S := □ 

xeM 

Let us assume that the hermitian structure is of the form H = rj + ioj, where rj is a 
Riemannian structure and oj is a complex structure. To define the measure and the volume 
form we can use either the complex structure oj (which defines a measure on T*M) or the 
Riemannian structure r] (which defines a volume form on the tangent bundle TM). 
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6.3 Average Operator Acting on Sections 

The averaging operation can be extended to a family of operators acting on sections of tensor 
bundles. This is especially important for the next section. In particular, let tt*, 7Ti, 7T2, 7r*T^ , ' 9 ^M 
and T(P'«)M as before. Let us consider the sections S G r(T(P'«)M) and ir*S G r^T^M) 
and the following isomorphisms: 

7r 2 | : T(ir* T (p,?) M ) — > r(T^M), S w i — ^ 5 Z , 

vr* : r(T( p '^M) — > r(7r*T^M), S 2 ^ <5 2 . 
Both isomorphisms are defined pointwise 

Definition 6.2 Consider the family of automorphims 

{A(W) : r(vr*TM) — »• r(VTM), W G 7r _1 (U), U G M}. 

The averaged operator of this family is the automorphism 

<A>: T(TuM) — »• r(TuM) 

suc/i i/iai at each point x G U it is given by: 

u G vr" 1 ^), S G TTM, 
where [Air* ■ S)(u) is the evaluation of the section [Air* ■ S)(u) at u. 

A similar definition application if the operators act on cartesian products of T (^T&'^M.) . 

7 Averaged Connection of a Linear Connection on 7r*TM 

We adopt a differential volume form uj x (y) such that {duo x {y))\s x = 0. Therefore we denote 
(wx(v))\v x = dvol(x,y). 

Definition 7.1 Let M be a n-dimensional smooth manifold, ir(u) = x and consider a differen- 
tiable real function f G JM. Then tt* f G .F(N) is defined by the condition 

<f = f(x). (7.1) 

The horizontal lift using a non-linear connection defined on the bundle TN — > N is 

i : TTM — ► TTN 

Xl ^ x ^ xl ih^ "^"^ (7 - 2) 

Proposition 7.2 Let M be a n-dimensional manifold and assume that N is endowed with a 
non-linear connection, u G ir~ l {x) C N, with x G M and let us consider a linear connection 
V defined on the vector bundle 7r*TM — > N. Then there is defined on M a linear covariant 
derivative along X, < V >x determined by the following conditions: 
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1. VX G T-jM and Y G TM, the covariant derivative ofY in the direction X, is given by 
the following averaging operations: 

<V> X Y :=< ir 2 \ u V Lu{x) ir* v Y > u , G U„, (7.3) 

where JJ U is a open neighborhood of u. 

2. For every smooth function f G J-~M the covariant derivative is given by the following 
average: 

< V > x f ■=< tt 2 | u V \ u{X )Kf >u, G U„. (7.4) 

3. The distribution function is such that < f > < oo. 

Proof: it is shown in reference [7, section 4}- a 

For the physical example that we are interested the averaging operation and the notion of 
volume that we are using is obtained by isometric embedding of the ambient Lorentzian structure 
r] on the unit hyperboloid times the positive weight function / (which will be the one-particle 
distribution function). Since in particular the function f x will be required later to be L 1 (S X ) 
and with compact support, 

vol(Y^ x ) := / f(x,y)dvol(x,y) < oo. 

The manifold ~E X is oriented. In particular, the integration is performed in the unit tangent 
hyperboloid, 

:= {y € T X M, | n(y,y) = l,y°> 0}. 

The measure that we will use is a solution of the Vlasov equation, which is obtained for a 
particular Liouville vector field. 

Definition 7.3 (Generalized Torsion) Let V be a linear connection on 7r*TM — > N ; then the 
generalized torsion tensor acting on the vector fields X, Y G TM is defined as 

Tor : r-7r*TM x r7r*TM — > IVTM 

(tt*X,tt*Y) -^Tor(V)(ir*X,TT*Y) = V l{x) tt*Y - V l(y) tt* X - tt*[X,Y}. (7.5) 
This tensor mimics the usual torsion tensor Tor, 

Tor : TTM x TTM — > TTM 

(X, Y) — > Tor(V)(X, Y) = V X Y - V Y X - [X, Y]. (7.6) 



Proposition 7.4 The averaged connection < V > has a torsion Tor(< V >) such that 

Tor(< V >) =< Tor(V) > . (7.7) 

Proof: It is shown in reference [7]. □ 

Corollary 7.5 Let M be an n- dimensional manifold and V a linear connection on the bundle 
tt*TM — > M with Tor(V) = 0. Then Tor(< V >) = 0. 
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Proof: It is direct from the proof of proposition (3.3.4). □ 
If Tor((< V >)) = we say that the connection (< V >) is torsion free. 

Corollary 7.6 Let M be an n- dimensional manifold. If the connection V in 7r*TM has the 

connection coefficients Then the average connection < V > has the coefficients 

< Fjk > (x) = 1 / T i jk (x,y)dvol(x,y). (7.8) 

V0L{2j x ) J-£ x 

Proof: Let {ei}, {7r*ej}, {i{ei)} be local frames for the sections of the vector bundles TM, 
7r*TM and the horizontal bundle % such that the covariant derivative is defined through the 
relations 

V^)7r*e fe = T l jk 7r*ei, i,j,k = 1,2,..., n. 
Then, let us take the covariant derivative 

< V > e e k = 1 v f 7T2(V t(e )n*e k ) = j [ ir 2 T) k ir* ei = 

1 



uoZ(S x ) Jy, x 

The relation (7.8) follows from the definition of the connection coefficients of the averaged 
connection, 

< V > e . e k =< r >* jk (x)ei(x) 

□ 



8 The Averaged Lorentz Connection 

Given the linear connection L V on the bundle 7r*TM — > S we can obtain an associated 
averaged connection. 

Usually the measure is given as f(x,y)dvol(x,y), where f(x,y) is the one-particle probability 
distribution function of a kinetic model. It must be gauge invariant and such that the low order 
momentum moments must be finite. The volume form is 

dvol(x,y) = v ^de~U ] ^dy 1 A • • -dy n ~\ y° = y°(x°, x n ~\ y\ .., y^ 1 ). 
Then one can prove the following: 

Proposition 8.1 The averaged connection of the Lorentz connection L V on the pull-back bundle 
7r*TM — > I] is an affine, symmetric connection on M. The connection coefficients are given 
by the formula: 

< L r jk >:= jk + (F l , < l -y m >r]mk + F\< X -y m > r, mj )+ 

+ F i m ±(< y m f 2 > r) jk - n^M < y m y s y l > ). (8.1) 
Each of the integrations is equal to the y -integration along the fiber: 

vol(S x )= j f(x,y)dvol(x,y), <y l >:= — -j— [ y % f(x,y)dvol(x,y), 
is, uoZ(S x ) Jv x 

< y m y s y l >:= 1 / y m y s y l f(x, y) dvol(x, y). 
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Proof: Equation (8.1) follows easily from the definition of the averaged connection by linearity 
We only need to prove that < y l > and the other moments are given by the corresponding 
integrals and the identity operator Id : 7r*TM — > 7r*TM, (x, y) >->■ y): 



<y l >= — 77=-: [ TT 2 Idy l TT*(y' l )f(x,y)dvol(x,y) = — -!— / y l f(x,y) dvol(x,y) 

VOL{2j x ) J-£ x V01(2j x ) J-£ x 



and similarly for other moments. Note that since y £ S x , rj(y,y) = 1 and it does not appear in 
the moments. □ 

Remarks. 

1. If we consider the bundle 7r*TM — > N, the coefficients of the averaged Lorentz connection 
are 

< L V >:= "Fkjix) -rf\x)(dA) lm {x) < ^-J-^ Vjs ( x ) y s y m Vkp (x) y? > + 

+rf l {x){dA) lm {x) < - 7= - Vjk(x)y m > +r, ll {x) < (dA) lk (x)— T ==-r jjs y s > + 

+rf\x){dA) lj {x) < -= -. m.s(x)y s > • 

2. The definition given by formula (8.1) does not depend on the natural coordinates. 

The following proposition enumerates basic properties of the averaged connection. The proof 
is straightforward; one needs to check the coefficients given by the formula (8.1): 

Proposition 8.2 Let (M, r/, [A]) be a semi-Randers space, f : XI — > R which is non-negative 
with compact support and < L V > be the averaged Lorentz connection. Then 

1. < L V > is an affine, symmetric connection on M. Therefore, for any point x 6 M, there 
is a normal coordinate system such that the averaged coefficients are zero. 

2. < L V > is determined by the first, second and third moments of the distribution function 
f(x,y). 

Remark While the first property is a general property of the averaged connection, the second one 
is a specific property of the averaged Lorentz connection and that we are considering trajectories 
whose velocity fields are in the unit hyperboloid X. 



9 Comparison Between the Geodesies of L V and < L V > 
9.1 Basic Geometry in the Space of Connections 

Let us consider the Lorentzian manifold (M, n) with signature (+,—,...,—) and M n-dimensional. 
If there is defined on M a time-like vector field U normalized such that rj(U, U) = 1, one can 
define the Riemannian metric fj [19]: 

fj(X, Y) := - V (X, Y) + 2 V (X, U) V (Y, U). (9.1) 
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r\ determines a Riemannian metric on the vector space T X M rj x = r]ij(x)dy % <8> dy J . Therefore 
the pair (T X M, fj x ) is a Riemannian manifold. It induces a distance function d^ on the manifold 
TyVI, df, : TyVl x T X M — »• R; (y,z) ^ inf{J a r)x{v,&)dT, a : I — >■ T,M, <r(0) = 
y, a(l) = z). 

We assume that /(x, y) has compact support on the unit hyperboloid bundle XI. Then the di- 
ameter of the distribution f x := f(x, •) : S x — > R; (x, y) ^ /(x, y) is a x := sup{d f f(y 1 ,y 2 ) \yi,V2 G 
supp(f x )}. We define the parameter a := snp{a x , x G M}. 

We choose as vector field [/ in the definition of the Riemannian metric (9.1) the following: 

U( X ) = { v»C)^£U>W i« »>.<»»(*)>» 

[ 0, otherwise. 

?7(x) is not continuous on the boundary d(^Tr(supp(f))), since in the interior of ir(supp(f)) is 
time-like and outside the closure of ir(supp(f)) is zero. In particular one has, 

H / x,| _ / rj(<y>,<y >)(x) > 1, x G int(supp(ir(f(x,y))) 

\\<y>(X)\\v I 0j x G M \ int(supp(ir(f(x, y))). 

One can approximate the main velocity field < y > by another vector field e < y > (x) which is 
smooth in the whole space-time M and time- like in an open sub-set of the interior of ir(supp(f)) 
(basically in n(supp(f)) \ dw(supp(f))). In this way, one can extend the definition of the vector 
field U (x) to the whole supp(f) in a smooth way. 

Proposition 9.1 Let the distribution function /(x, y) be of class C 1 and rj a C°° Lorentzian 
metric on M. Then there is a C 1 vector field 6 < y > which coincides with the time-like vector 
field < y >, and such the following relation holds: 

\\ <y> - <y>\\e fj = 0{a 3 ), 

with a the diameter of the support of the distribution f x and with e fj the associated Riemannian 
metric to the vector field 6 < y >. 



Proof: Let us consider the C°° function 

b(t) = 



0, t < 



k is a positive parameter. Now consider an arbitrary open subset e O C int(supp(ir(f(x,y))) 
and the restriction of the vector field < y > to e O denoted by < y > \ e£>- Then, we consider 
the distance function between each point z G d( y supp7r(f)) to the open set e O. This distance 
function, which we denote by t(z) is realized by a point w G e O. We can consider the geodesic 
segment between z and w, for a similar enough e O to int(supp(Tr(f(x,y))) in such a way that 
these segments exists. After this, we can define the following interpolating vector field: 



c -v(x) 



0, x G M \ int(supp(7r(f(x, y))), 

t(z)v(w(z)), x = z G int(supp(ir( f(x, y))) \ e {O}, 
v(x), x G e O. 



This vector field is still not smooth but it is continuous. One can smooth it using bump 
functions based on the function b(t) in such a way that the resulting smooth field is zero in 
M \ int{supp{T\-{f{x,y))) and where the parameter t depends on the point z G d( y suppir{f)^ 
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and essentially corresponds to the geodesic distance (using 77, which we now it is smooth on 
int(supp(iv(f(x,y))). 

The construction of the metric e f] is based in the fc-parameter family of bump functions based 
on the function b(t). The idea is to fit the parameter k such that in the smoothing process using 
bump function [27], we obtain || < y > — < y > \\ Cf j = 0(a 3 ). The construction is similar to 
Steenrod's construction [39, section 6.7]. □ 

Given a continuous operator A x : T X M — > T^M, its operator norm is defined by 
\\AUx) := sup { yeT x M\ {0} }. 

\\y\\r] J 

Let us denote the space of linear connections on TS by Vs. This space is a finite dimensional 
manifold which points are coordinated by the set of functions {P^(x, y), k = 0, 1, ...,n — 1}. 



Proposition 9.2 On the space Vs, there is a distance function. The distance between two 
points 1 V, 2 V G Vs is given by: 



1 V, 2 V G V s }. (9.2) 

Proof: The function (9.2) is symmetric and non-negative. The distance between two arbitrary 
connections is zero iff 

v / r?( l V x X - 2 V X X, ^V X X - 2 V X X) = 

for all X G ~E X . This happens iff ^ xX = 2 VxX for any X G S x . The triangle inequality also 
holds, using the triangle inequality for 77: 



d,(V 1 ,V 3 ) = su P { }< 

< sup I ■ 
f 





2 VxI, X V X ^ - 


- 2 V X X) 


Vv(x,x) 


0?( 2 V X X- 


3 V X X, 2 V X I - 


3 ViI) ] 



^d^V, 2 V) + ^( 2 V, 3 V). 



□ 



For an arbitrary 1-form ui we denote by oj^ := r/ -1 ^, •) the vector obtained by duality, using 
the Lorentzian metric n. Similarly, given a vector field X over M, one can define the dual one 
form X* := n(X, •); lxlo is the inner product of the vector X with the form u. 

Proposition 9.3 Let f(x,y) be the one-particle probability distribution function such that each 
function f x has compact and connected support supp{f x ) C T, x . Then 

( L V y y- < L V > y y)(x) = -(i S F)t(x) • (t y (S)(x,y)) + 2 (5 2 (y))(x,y) + 3 (5 3 (y))(x,y), (9.3) 
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where S m (x,y) =< y m > (x) — y m does not depend on y £ supp(f x ) or the 2-form F. The 
tensors Oi(x,y) are given by the following expressions: 

0?(5 2 (y))(x,y) = ^<y m > ( x )S s (x,y)S l (x,y)+ < y m X 5 s (x,y)6 l (x,y) > 

+ 2<y l >(x)< 5 m (x,y)5 s (x,y) > )r ]am ^y k , (9.4) 
Of{5\ y ))(x,y) = \ < 5 m (x,y)5 s (x,y)5 l (x,y) > r, 8jmk y> i y k '. (9.5) 
Proof: From the expressions for the connection coefficients, 

< L V > y y - L V y y = l(F^(x)(< y m > (x) - y m )n mk + F* fc (x)(< y m > (*) - y m )Vm 3 )+ 

+F* m (x){(< y m > (x) - y m ) Vjk - n ]S r, kl (< y m y s y l > (x) - y m y s y l )))v V 

since rnjy % yi = 1. The difference between the two connections can be expressed in terms of the 
following tensors: 

5 m (x,y) =< y m > (x) - y m , 5 msl (x,y)y syi =< y m fy l > (x) VsjVlk y> y k - y m 
and is given by the following expression: 

(< L V > y y- L V y y)(x) = ^(F l 1 (x)(5 m (x,y)) Vmk +F l k (x)(5 m (x,y)) Vmj + 

+ F i m (x)(S m (x,y)r ljk -r ljsm d^(x,y)))y^y k -^- = 

= (F* jy i6 k (x, y)y k + F* m (x) (5 m (x, y) - y 3 y k 5 mk \x, y)))^, 
where yj = n jk y k . 

In the above subtractions the second contribution is of the same order in S(y) as the first one. 
To show this, recall from the definitions 

5 ms \x,y)y syi =< y m (x)y%x)y l (x) > W yV - y m 

where the hat- notation is used to distinguish integrated variables from the coordinates y % . Then 
we can use the following relations: 

f =< f > (x) - 5(x,y s ), f =< y s > (x) - 5 s {x,y), y m =< y m > (x) - 5 m {x, y). 

One substitutes this relation and, taking into account that < S(x, y) >= 0, one gets 

5 msl (x, y)y s yi = ( < y m > (x) < y s > (x) < y l > (x) + < y m > (x) < 5 s (x, y)5 l (x, y) > + 

+ < y l > (x) < 5 s (x,y)5 m (x,y) > + < y s > (x) < 5 m (x,y)5 l (x,y) > - 
- < 5 m (x,y)S s (x,y)5 l (x,y) > ^jr] sj r]i k y j y k - y m . 
Now we use a similar relation to go further: 

y s =< y s > (x) - 5{x, y s ), y l =< y l > (x) - 5 l {x, y). 
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We introduce these expressions in the calculation of 5 m& : 

5 msl (x,y)y syi =(<y m >(x) (y s + 5 l (x,y))(y s + S l (x,y))+ < y m > (x) < 5 s (x,y)5 l (x,y) > + 

+(y l + 5 l (x,y))<5 s (x,y)6 m (x,y)> +(y s + 5 s (x, y)) < 5 m (x, y)6 l (x, y) > - 

- < 5 m (x,y)S s (x,y)S l (x,y) > ^r] sj 7] lk y j y k -y m . 

Using again the fact that y l r\\ k y k = 1 we get (again using that < y m >= y m +S(x, y) to recombine 
the first and last term): 

S msl (x,y)y syi = 5 m (x,y) + 2 < y m > (x) 5 s (y)^ + 20?(5 2 ) + 20 3 m (<5 3 ). 

The tensors Oi are given by the formulas (9.5) and (9.6). 

The transversal contribution to the difference between the connections is given by: 

\^ m {x)(5 m {x,y) V3k - msm 5 ms \x,y)yiy k ) = l -F* m (x) (5 m (x, y) -5 m (x,y) 

-2<y m > (x)6 s (x,y)r) sj yi) - 0?(5 2 ) - 3 m (<5 3 ) = 
= -F < m (ar)( <y m > {x)5 s {x,y)n s ^) - 0^{5 2 ) - 3 m (5 3 ). 
The longitudinal contribution is 

i(F* j(x)(< y m > (x) - y m ) Vmk + F\(x)(< y m > (x) - y m )v mj ) V 3 y k = F* j(x) y'd k (x,y)y k . 

Adding together the longitudinal and transversal contributions and taking into account the 
formula 5 m (x,y) =< y m > (x) — y m we get the following expression: 

(F i jy ^)(S k (x, y )y k )-(F i m ( X ){<y m >(x)d s (x, y ) Vsj ^) + 2 {5 2 ) + 3 (5 3 )) = 

= -(F l /(^)) (5 k (x,y)y k ) - 0%{5 2 ) - 3 m (<5 3 ). 

□ Let us consider a frame {ej, i = 0, n — 1} such that eta is diagonal at the point x G M. 
After calculating the distance of the connections using the formula (4.6.3), the leading term in 
5 is quadratic. 

Proposition 9.4 Let (M, rj, [A]) and L V be as before and assume that f x has compact and 
connected support for each fixed x G M and that a := sup{a x , x G M} << 1. Then the 
following holds: 

(%( L V, L V >){x) < ||F||^(x)C(x)q 2 + 2C 2 2 {x)a 2 {l + a) + C 3 3 (x)a 3 (l + a), (9.6) 
with C(x), C2(x), C 3 (x) being functions depending only on x with value of the order of unity. 

Proof: From equation (4.6.3) one obtains: 

|| L V y y- < L V > y y\\ f} = \\F i j(x)P(x,y)5 k (x,y)y k e t +0 2 (5 2 )(x,y) + 3 (S 3 )(x, y)\\ r , < 

< \\F i jP(x,y)5 k (x,y)y k e l \\ n +\\0 2 (5 2 )(x,y)\\ n + \\0 3 (5 3 )(x, y)\\ n . 
Each of these three terms can be bounded. 

The bound of the first term comes as follow, since ||ej|| = 1: 

\\F i j5 j (x,y)5 k (x,y)y k e i \\ fl < ||F||^ • \5 k (x,y)y k \ < ||F||^ • \\S(x, y )\\fj ■ \5 k (x,y)y k \. 
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Since the support of the distribution function f(x,y) is compact and connected, one can write 
the decomposition < y > (x) = e(x) + z(x) with the property that z(x) G supp(f x ). In the case 
of the Minkowski metric one can check by geometric inspection that ||e(x)||»? < a. 

This bound of e(x) follows from the shape of the unit hyperboloid in the following way. First, 
note that the domain S x := {y £ T^M | r/(y, y) > 1 } is a convex set. Indeed, we note that 
d(J} x ^ = H x and that < y > (x) 6 S x . Secondly, each (Yl x ,f] x ) is a Riemannian structure with 
zero curvature tensor. Therefore it is isometric to a piece of the Euclidean space. Therefore, we 
can use the standard definition of center of mass, in this case with an invariant measure given 
by f(x, y, s) dvol(x, y). The function f(x, y, s) is such that 

[ dsf{x,y,s) = f(x,y), 
J o 

where s is the parameter of the convex line connecting y £ supp(f x ) with < y >. Let us denote 
by swpp(f x ) the convex envelope of supp(f x ). By construction < y >G supp(f x ). Indeed one 
can check that < y > is the center of mass of the convex set swpp{f x ) [36]. 

We have the following bound 

\\d(x,y)\\fj < || < y > (x) - y\\fj < \\e + z(x) - y\\ n < 

^ \\ e \\v + \\ z ( x ) ~ y\\fj < a + a = 2a. 
For the third factor, one has the following bound (5(x,y) =< y > —y): 

\S k (x,y)y k \ = \<y k > (x)y k - 1| = | < y k > (x)(y k - <f>(x)+<y k > (x)) - 1| < 

< | < y k > (x)(y k - <y k > {x))\ + | < y k > (x) < y k > (x) - 1|. 
Using Cauchy-Schwarz inequality for f) we obtain: 

\S k {x,y)y k \ < || < y k > (x)\\ n \\{y k - < y k > (x))||^ + | < y k > (x) < y k > (x) - 1| < 
< || < y k > (x)\\ a + | < y k > (x) < y k > (x) - 1| < Jl + j^a + (Jl + ||e||^ - 1) < 



< VTT~aa + + a - 1) < (1 + a)a + (1 + a - 1) = 2a + a 2 . 
Therefore, one obtains: 

\\{^^(x,y)){5 k (x,y)y k )\\^ < \\FUx)C(x)a 2 . 

The function C(x) in equation (4.6.6) is bounded by the constant 2. 

Using homogeneity properties, one can see that the following relations hold: 

\\02(b 2 )h < C 2 2 (x)a 2 {l + B 2 {x) a), (9.7) 

and 

||0 3 (5 3 )|^ < Cl(x) a 3 (l + B 3 (x) a). (9.8) 

The functions Cj(x) depend on the particular shape of the support of the distribution function 
/ and on the curvature of the metric r\. Using geometric arguments (and in particular compact- 
ness and connectedness of the supp(fx)), one can bound the function \S°(x,y)\ in terms of a in 
a similar way as we did for C(x). The reason why the constants will be of order 1 is that this 
was the case for C{x) and there is not new divergences in the bound of the functions Bi(x) and 
C ( x). □ 
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Remark 1. For arbitrary Lorentzian manifolds the same ideas can also be applied to bound 
the functions Cj(x), as we have directly seen for C{x). A homogeneity argument implies that 
these functions are also bounded by constants of order 1. 

Remark 2. In this proof it is not essential that a « 1. However it simplifies the calculation. 

Corollary 9.5 Let (M, i], [A]) be a Randers space with U^eM su PP{fx) compact, with L V and 
< L V > as before. Then there is a global bound: 

d n { L V,< L V>)(x)< C||F||^« 2 + 2C 2 a 2 (l + B 2 a) + Cfa 3 (l + B 3 a), Vx G M. 

where the constants C, C 2 , C3 , B 2 , B$ are of order 1 . 

Proof: It follows from proposition (9.4) and compactness of IJzeM su PP(fx)- D 
9.2 Comparison Between the Geodesies of L V and < L V > 

There are several ways of defining the energy of a bunch of particles. We have choose one which 
will be useful to get the desired comparison. We define the energy function E of a distribution 
/ to be the real function: 

E : M — ► R 

x ' V E(x) := inf{y°, y G supp(f x )}, (9.9) 

where y° is the 0-component of a tangent vector of a possible trajectory of a point particle, 
measured in the laboratory coordinate frame. 

Let us restrict ourselves to the case that the Lorentzian metric is the Minkowski metric in 
dimension n. Let us denote by j(t) the gamma factor of the Lorentz transformation from the 
local frame, defined by the vector field U to the laboratory frame, at the instant t defined by the 
local time, the coordinate time defined by the laboratory frame. Let us define 2 (t) = y 2 ^)— < 
y > 2 (t) and 2 {t) =< y > 2 (t) — y(t). Here y(t) is the velocity tangent vector field along a 
solution of the Lorentz force equation and y(t) is the spatial component of the tangent vector 
field along a solution of the averaged equation, with both solutions having the same initial 
conditions. The maximal values of these quantities on the compact space-time manifold are 
denoted by 8 2 and 9 2 . 

Theorem 9.6 Let (M, 77, LA]) be a semi-Randers space and rj the Minkowski metric. Let us 
assume that: 

1. The auto-parallel curves of unit velocity of the connections L V and < L V > are defined 
for time t, the time coordinate measured in the laboratory frame. 

2. The ultra-relativistic limit holds: E(x) » 1 for all x € M. 

3. The distribution function is narrow in the sense that a « 1 for all x G M. 
4- The following inequality holds: 

\6 2 - 9 2 \ < 1, 

5. The support of the distribution function f is invariant under the flow of the Lorentz force 
equation. 

6. The change in the energy function is adiabatic: j^logE « 1. 
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Then for the same arbitrary initial condition (x(0),x(0)), the solutions of the equations 

L V ± ± = 0, < L V >~ x = 

are such that: 

\\x(t)- x(t)\\ < 2(C(x)||F||(x) +C|(x)(l + B 2 (x)a))a 2 E~ 2 (x)t 2 , (9.10) 

where the functions C(x), Ci(x) and Bi(x) are bounded by constants of order 1 and the norm 
|| • || is the spatial norm in the lab. frame. 

Proof: At the instant t we calculate the distance measured in the laboratory frame between x(t) 
and x(t), solutions of the geodesic equations of the corresponding connections L V and < L V >, 
when both geodesies have the same initial conditions (x(0),x(0)). For this, we use the general 
formula for the solution: 

x i (t) = x i {d) + J^ds(x{d) + dlx\l)y (9.11) 

Since the initial conditions for both connections are the same, the equivalent relation for the 
geodesies of the averaged connection is 

x l {t) = x*(0) + J ds(x(p) + J dW{lj} . (9.12) 

We have to estimate the distance between both solutions. Since we know the distance between 
both connections, it is possible to give a natural bound for the distance between solutions. The 
main property is the smoothness theorem on the dependence of solutions of differential equations 
on the external parameters (for instance [20, Appendix 1] or [30, chapter 1]). 

Let us consider the family of connections depending on the distance between the two connec- 
tions 

Uax = d n ( L V,< L V >) 

given by the convex sum: 

f V := -^(W - £) L V +j^—Z< L V>, £€[(Umax]- (9.13) 

srnax srnax 

For £ = one has = L V, while for £ = ^ m ax one has the averaged connection. Using 
the result of the smoothness of the solutions of the differential equations, one can expand the 
solution ^x l of the geodesic equation for the connection with parameter £ in the following way. 
In particular, the second derivative with respect to the coordinate time t reads: 

= (d^x% =0 ^ + O(e). (9.14) 

We need to evaluate the derivative (d^x l )\^ = o. From the formula (4.6.14), one obtains: 

(d^x%=o = ^— ■ ( L V x x- < L V x x >)', 

smai 

where x(t) is the solution of L V±x = with the given initial conditions. This is because the 
support of the distribution function / is invariant under the flow of the Lorentz force equation. 
Therefore, since we are dividing by the distance £ max , the derivative is such that its norm is 1: 

||(9 C c £)|{=olh = T~ ■ IK L ViX- < L V ± x >)||^ = 1, 
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because of the definition of ^ mas an d the formula (4.6.14). On the other hand the relations 
between proper times and coordinate time in the laboratory frame are 

dr = 7~ 1 di, df = 7 _1 c/t. 

This implies the following relation between derivatives: 

d _ 1 d d d 
dt ^ dr' dt ^ df 

Using the hypotheses \6 2 {t) -9 2 {t)\ « 1 and j^logE « 1, one obtains the following relation, 

The factor 2 comes from the bound of the term where appears the derivative of the energy. We 
have use the adibatic hypothesis. 

Therefore, in the ultra-relativistic regime, 

\\x{t)-<t)h< 2a 2 E- 2 \\F\\ n \\-^— ■ ( L V Z Z- < L V z Z>)\\ n t 2 . 

£,rnax 

In order to transform the distance function based on the norm || • ||^ to the distance function 
based on the spatial norm || • || defined in the laboratory frame, one uses the Lorentz contraction 
factor 7, 

\\x{t)- X {t)\\=T 1 \\m-<t)u 

we obtain: 

\\x(t)- x{t)\\ < 2(C\\F\\ fj (x) +C|(x)(l + B 2 {x)a))a 2 ^ 1 (x)E~ 2 (x)t 2 . 

The norm of the electromagnetic field also changes when it transforms from one coordinate 
system to another coordinate system. The way of changing is given by a Lorentz transformation, 
which gamma factor is j. Therefore: 

l|F|| =7l|F|k 

Finally we get, 

\\x(t) - x(t)\\ < 2(C||F||(ar) + C 2 (x)(l + B 2 (x) a))a 2 E~ 2 (x)t 2 . 

□ 

Remark. It is important to notice that the above bound is Lorentz covariant. Although we 
have performed the calculation in a particular local coordinate frame, both the left and right 
sides of (9.10) obey the same transformation law under Lorentz transformations. For instance, 
if we express this bound in a coordinate system in the co-moving frame of the bunch, one has 
that: 

\\x(t)- x(t)\\a< 2(C||F|| j j(a;) +Cl{x){l + B 2 {x)a))a 2 m~ 2 P, 

for a mass m = 1 in the units we use. One can check that formula (9.11) is Lorentz invariant 
by direct inspection of the transformation rule of each of the factors that appear there: 

1. o is Lorentz invariant. 

2. The product E~ 2 t 2 is Lorentz invariant. 
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3. The functions C(x), C\(x) and B\{x) are Lorentz invariant. 

4. The function ||F||^ transforms in the same way as \\x(t) — x^)^ 1 under Lorentz transfor- 
mations. 

5. The condition dlo f t E << 1 does not appear explicitly in expression (9.10). Therefore we 
do not need to check it explicitly. 

Corollary 9.7 Let (M, 77) be as before and such that there is a global bound for ||F||^(x) < 
||F||f7 < oo, the energy function is bounded from below by a constant E and the curves x(t) and 
x(t) are compact. Then there are some constants Q such that Ci(x) < Ci and the following 
relation holds 

\\x(t)- x(t)\\ < 2(C(x)\\F\\(x) +C|(1 + a))a 2 E~ 2 t 2 . (9.15) 

Remarks 

1. In the above result the 2-form F is physically interpreted as the Faraday form. 

2. These bounds happen at least in two physical scenarios: 

(a) That the manifold M is compact. In this case, one has to consider spatial boundaries, 
since it is well-known that if the space-time is compact, they exists closed time-like 
curve ([34, pg 58]). This can broke up causality. 

(b) The trajectories that we consider are compact. In this case, one can define an ef- 
fective compact space-time manifold with a spurious boundary and apply the first 
case, although excluding closed time-like curves, which do not happen in physical 
situations. 

3. The external field F does not depend on the energy of the beam of particles. However, this 
is not necessarily the case in some situations. For instance, there are effects which could 
be reduce the beam size (adiabatic damping and Landau Damping [10]). If this happens, 
there is a strong reduction of the size in energy and momenta of the beam. This implies 
that one can describe this as an effective exponent E~ 2+ ^ , with j3 negative. 

4. That the curves x(t) and x(t) have compact image has a physical interpretation: all the 
trajectories start in a source region and finish in a target region. 

Theorem 9.8 Under the same hypothesis as in theorem 4-6.6, the difference between the tangent 
vectors is given by 

\\x(t) - x(t)\\ < (K(x)\\F\\(x) + Ki(l + D 2 (x)a))a 2 ) E- 1 1, (9.16) 
with Ki and D2(x) functions bounded by constants of order 1. 

Proof: The proof of this theorem is similar to the proof of Theorem (9.6), although based on 
the following formula for the tangent velocity field along a curve: 

x(t) = x{0) + f x(l)dl. (9.17) 
J 

□ 

By the same argument the relation (4.6.16) is Lorentz invariant. 

Corollary 9.9 Under the same hypothesis than in Corollary 9.7, there are some constants of 
order 1, K , Ki and D, L , such that K(x) < K, K(x)i < Ki, Di(x) < Di and the following relation 
holds: 

\\x(t) -x(t)\\ < (K\\F\\(x) +Kl{l + D 2 {x)a))a 2 )E- 1 t, (9.18) 
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10 Applications of the Averaged Lorentz Dynamics and Discus- 
sion 

Theorem 9.6 and Theorem 9.8 are proved under some hypotheses which have to be motivated 
from a mathematical or physical point of view. 

1. The first hypothesis asserts the existence of solutions of some geodesic equations at least 
in an interval of time [0, T] 3 T. If the space-time manifold M is compact and with 
boundary, this hypothesis is strong, since there is a time T where any geodesies finish in 
a boundary. If we are in a situation such that the solutions of the auto-parallel equations 
are compact sub-sets on M, we must guarantees the existence of the solution up to time 
t < T. 

2. Ultra-relativistic dynamics implies (as an assumption) that the minimal value of y° on the 
support of the distribution is much larger as 1. The definition of ultra-relativistic limit is 
a gauge invariant notion and it is also independent of the frame of reference. 

3. The condition that states that \6 2 (t) — 6 2 (t)\ <C 1. If the motion is ultra-relativistic this 
condition implies = 1 + e, with e << 1. Therefore the above hypothesis is the 
statement that the integral curves corresponding to the averaged connection are in the 
ultra-relativistic regime too. 

4. A similar interpretation for the narrowness condition is clear: the diameter of the distri- 
bution is much smaller than the energy of the system. The condition of narrowness is 
gauge and Lorentz invariant. We can compare our definition of narrowness (which is that 
E{x) » a) with the invariant warm fluid condition, which is that [29] 

< (y k - < y k >){y j - < y j >) > Vjk « 1. 

That condition is Lorentz invariant and provides a notion of narrowness, since the quantity 
< (y k — < y k >)($— < y-* >) > is just the variance of the distribution /. 

5. Our condition \9 2 {t) — 2 (t)\ <C 1 is equivalent to the invariant warm fluid condition, 
because < (y k — < y k >)($— < y-' >) >= \9 2 {t) — 2 {t)\ when everything is applied on 
the unit hyperboloid. 

6. The hypothesis on the invariance of supp(f) by the flow of the Lorentz force equation 
guarantees that the distribution functions that we consider are related with the Lorentz 
force equation and therefore that they are physically significant. One possible relation is 
through the Vlasov equation, which is an equation describing a distribution function in 
kinetic theory [18]. In Vlasov's model, the one-particle probability distribution function f is 
constant along the flow of the Lorentz force. From this perspective, this is the most simple 
case, since it could be that supp(f) is invariant by the flow but not constant. Therefore, 
our theorems (9.6) and (9.8) apply to more general kinetic models. 

There are two situations where the results of section 9 are applicable: 

1. The ultra-relativistic regime, E — > oo. In this case, \\x(t) — x(t)\\ — > for distributions 
with finite a and for finite times t. The requirement of \9 2 — 9 2 \ ^ 1 is naturally fulfilled 
for common distributions and in particular for smoothly truncated Gaussian distributions. 
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2. For the Dirac delta distribution: 



f(x,y) = 5(y-V(x))V(x), x G M. 



(10.1) 



This delta distribution function corresponds to the charged cold fluid model solution of 
the Vlasov equation. 

Since the width in the space of tangent velocities of the distribution is zero, a = 0. 
The trajectories, however, are constrained to follow the Lorentz force equation, since by 
hypothesis the support of the distribution function is invariant under the flow of the Lorentz 
force equation. Note that in this example, the distribution is not smooth. However, one 
can check that the Lorentz and averaged Lorentz connections coincide. This fact suggests 
that we can extend our results to a bigger category of functional spaces for the distribution 
function /. 

Let us discuss the limits of applicability of the averaged dynamics and in particular of the 
formula (9.10). We make the assumption that during the time of evolution the 7 factors are 
constant. This is not realistic in general (since it can happens acceleration of the beam of 
particles) but it is useful in order to provide an estimate of the range of validity of the time t of 
the approximating the Lorentz dynamics by the averaged dynamics. Therefore, equation (9.11) 
takes the form 



Let us assume a natural maximal distance Lq between two particles in a bunch. For instance, in 
an accelerator machine Lq could be the diameter of the pipe. The averaged model is applicable 
until the difference ||£ l (i) — x*(t)|| is of order Lq. This can start to happen after a time evolution 
such that the difference \\x l (t) — x l {t)\\ = Lq. The characteristic time where the averaged model 
loses validity is: 



Corollary 10.1 In the limit E — > 00, -^logE << 1 and (or) the limit a — > if all the other 
parameters remain finite, t max — > 00 . 

Applications of the Averaged Lorentz Dynamics 

There are some applications of the averaged Lorentz force equation: 

1. The solutions of the averaged dynamics can be used as a definition of the reference trajec- 
tory in beam dynamics. There are two reasons for adopting this perspective: 

(a) As Theorem 9.6 states, for the same initial conditions, in the ultra-relativistic limit 
and for narrow probability distributions, the difference between the original trajectory 
and the averaged trajectory is small. 

(b) For narrow distributions in the ultra-relativistic limit, the integral curves of the av- 
eraged velocity field coincide at least up to order a 2 with auto-parallel vector field of 
the averaged connection. This is proved in the companion paper [9]. 

The reference trajectory is a fundamental notion of beam dynamics [10]. It is usually 
understood that the reference trajectory does not necessarily coincide with the trajectory 
of a real particle. Our criticism of this idea comes from the point that in this case, the 
reference trajectory is not observable. Therefore, referring all the optic design and other 



x(t)-x(t)\\ < Ca 2 E- 2 {to)\\F\\t 2 . 




(10.2) 
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parameters of an accelerator machine to this trajectory can be a dangerous strategy in 
extreme high density and energy conditions, since one is assuming an ideal behavior which 
in principle can or can not be the right description of the dynamics in such situations. 

Based on this, we propose that the reference trajectory is directly linked with a solution 
of the averaged Lorentz dynamics. This method has the advantage that because the 
solutions of the averaged dynamics are observable, its solutions can be directly tested and 
experimentally controlled. 

2. The averaged dynamics is a fundamental tool in a new justification of the use of fluid 
models as an approximation of a kinetic model. There are two reasons for this: 

(a) It is a general belief in beam dynamics and plasma physics that fluid models appear 
from kinetic theory as an averaged description (for instance a standard reference is 
[11], [29] and [31]). However, derivations are based on assumptions on the higher 
order moments of the distribution function. These are hypotheses which can or can 
not hold at very extreme conditions. We have approach that problem without making 
hypothesis on the higher moments. 

(b) Since the averaged Lorentz connection is affine and symmetric, one can make use of 
normal coordinates to simplify some calculations. This technical point is fundamental 
in our approach to the problem. Indeed, this is a consequence of the fact that an 
averaged dynamics is technically simple compared with the original dynamics. 

3. In order to write down the geodesic equations of the averaged dynamics, one only needs 
the first, second and third order moments. This makes the averaged dynamics insensitive 
to higher order moments of the original distribution /. We think this is an advance which 
respect to standard treatments [11], [29], [31]. The approximation of the Lorentz dynamics 
by the averaged dynamics requires narrowness of the distribution function and that the 
process happens in the ultra-relativistic regime. Since the averaged model is insensitive 
to higher moments, it is not harmful to assume that moments of higher order as three are 
zero, as soon as one proves this is consistent with the premises of the approximation. We 
think that this is a general principle by which to understand the validity of fluids models 
as an approximate description of kinetic models. 

4. The averaged dynamics defined by the averaged connection < L V > can be used in 
numerical simulations of the bunch dynamics in accelerator physics. Since the averaged 
equation is simpler, it can be used to improve the calculational speed in simulations of the 
evolution of bunches composed by a large number of charged particles. This is an idea, 
however, that still needs to be tested. 



Relation of this work to previous work. 

The averaging procedure used in this work was first introduced in reference [7] by the author, 
although in the context of Finsler geometry and formulated using a geometric formalism. In 
this paper we have adopted a local description using frames and local coordinates. 

At first glance, there is a resemblance between the averaged procedure presented here and 
the procedure discovered by Zalaletdinov and co-workers [21]. However we think that these 
procedures are intrinsically different, although unfortunately the use of a common notation and 
terminology makes it easy to confuse them. In fact, the averaging considered by those authors 
is very different from the one here and also the one appearing in [7] ; their averaging procedure is 
taken on open sets of the space-time itself, while our averaging procedure is an integration along 
the fiber over a pull-back bundle over the unit hyperboloid. In local coordinates, it corresponds 
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to integration on the space of allowed velocities, at fixed point of the space-time. Indeed, our 
procedure has some resemblance to the averaging along the fiber in classical mechanics [22] or 
the averaging along the fiber that appears in algebraic topology, through a generalization of 
Poincare's lemma [23]. These two examples have been formalized in a general framework [8]. 

In this work we are considering a geometrization of the dynamics of charged point particles 
with external fields and not a geometrization of the whole electromagnetic theory. Indeed, one 
of the main points discussed in this work is that the natural interpretation of the classical 
electrodynamics of point charged point particles in an external electromagnetic field does not 
correspond to a space-time metric structure. There are theories which are a geometrization of 
Classical Electromagnetism and where Finsler and Lagrangian structures appear in a natural 
way. For such theories, one can see, for example, the approach due to R. Miron and co-workers 
[16], [24]. However, the aim of the present work is very different and must not be confuse with 
the kind of theories developed by these authors. 

There exists a previous geometrization of Lorentz force equation discovered by Miron and 
co-workers [37]. The connection that they define coincides with our V connection. We have 
seen that both are affine equivalent connections (ones the parameter has been fixed to be the 
proper-time parameter associated with f). However, we have not yet clarified what is the relation 
between the two approaches. 



11 Notion of Almost-Connection 

During our analysis of the Lorentz force equation, we have considered the associated non-linear 
connection on TN and the associated linear connections on 7r*TM — > S. However, strictly 
speaking the system of differential equations (1.1) does not define a non-linear connection on 
TS. The reason is the following. Let us fix a point u G S. There exists a natural embedding 
e : S ^-t- TM. Then one can consider the sub-bundle e(X!) and the covariant derivative L D 
acting on elements of rTe(Xl). This covariant derivative can be extended to derive sections of 
the extended bundle |_|ue£ T U N, with is a sub-bundle of TM: 

L D : r( | | T U N) x r( [J T U N) — > T( |J T U N), 

«es «es «es 

From the definition of L D, this operator cannot be extended in a smooth and natural way to 
be a covariant derivative on TM. This is because of the appearance of factors rj(y, y) in the 
connection coefficients: the function y/rj(y, y) is not defined in an open neighborhood. The 
non-extendibility of this function is the problem to extend the covariant derivative along those 
directions. 

The above fact suggests the existence of a mathematical object which is a generalization of the 
ordinary notion of covariant derivative in the sense that allows covariant derivatives along outer 
directions to the manifold, but that are not obtained as a restriction of an ambient covariant 
derivative. 

Example. The Lorentz connection was obtained from the Lorentz force equation assuming 
that the torsion is zero. We obtained the following connection coefficients: 

jk 1 jk ~r 1 jk ~r ■ lj jki 
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v m 1 \ 

T jh = F' m — — (7/j fc - — — 7VjsVkiy s y 1 ))- 

This defines a rule to derive sections of TE along directions of TS. The same coefficients 
provide a rule to derive sections of |J Mg £ T U TM along directions of TM, but we can extend the 
definition of the covariant derivative to the whole manifold TM. 

A related issue is the following. On the unit hyperboloid recall that T l jky^y k = (this is 
why it was called transverse component). Therefore, if one considers instead an alternative 
connection L V, determined by the connection coefficients 

L-pi r?-pi I rt 

1 jk — 1 jk "T J-i jk, 

the corresponding geodesic equation (parameterized by the proper time of the Lorentzian metric 
77) is again the Lorentz force equation [37]. Therefore we see that both V and V reproduce 
the Lorentz force and both are torsion-free (the connection coefficients are symmetric in the 
lower indices). This is in contradiction with the fact that a connection of Berwald type (linear 
or non-linear) is determined by the set of all geodesies as parameterized curves on the base 
manifold M and the torsion tensor (for linear connections the procedure can be seen in [17]; for 
non-linear connections, a procedure to define the connection is described for example in [16]. 

A solution to this dilemma comes from the fact that the Lorentz force equation applies only 
to time-like trajectories for which the tangent velocity vector fields live on the unit hyperboloid 
S; that is, we cannot extend the operator 

L D : r( [J TJM) x r( [J TJM) — > r( [J TJM) 

«G£ weSigma «ex 

to a genuine operator of the form 

L D : r(TM) x r(TM) — ► r(TM). 

In other words, we cannot extract enough information from the Lorentz force equation to de- 
termine a non- linear projective connection ([38]), because on the base manifold, we do not have 
information about all the possible geodesies. This is why we can not extend the connection. 

One can be tempted to extend the Lorentz force equation to another equation valid for any 
kind of trajectory. This is partially accomplished by the averaged connection. However, these 
kind of extensions could be non-natural or unique and indeed hide a natural object such as the 
almost- connection. 

Preliminary definition of Almost- Connect ion 

Let M be a manifold of dimension n. An almost-connection is a second order vector field 
X defined on a sub-bundle D TTM, with D of arbitrary rank. Associated with D is the 
corresponding almost projective covariant derivative on 7r*TM. 

Examples. 

1. A projective connection (in the sense of Cartan [38]) is an almost connection such that 
D = TM. 

2. The Lorentz connection L V provides an example where S = D / TM. One can define a 
Koszul connection acting on T( |J«es T U TM). 
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3. The averaged covariant derivative < L V > is an affine connection and therefore an almost 
connection in the above sense. 

One can consider the corresponding linear connections L V and L V on 7r*TM. Generally, their 
averaged connections < L V > and < L V > are not the same. However, we would like to have 
an averaged operation which is well defined for the objects in a given category. By definition 
this will be the category of almost-connections and the corresponding morphisms. We require 
that the result of the averaging operation be the same for each representative belonging to the 
same almost connection, 

Definition 11.1 Let M be a manifold of dimension n. An almost non-linear connection is the 
maximal set of semi-sprays \ defined on a sub-bundle D C TTM such that they have the same 
averaged linear covariant derivative < V > and the same torsion tensor: 

f(X,Y):= L b jt Y- L b f X-[X,Y], X,YeFT>. (11.1) 

Associated with x over D is the corresponding almost- covariant derivative on the bundle 7r*TM — 
D. Any of the connections in the same almost-connection has the same averaged connection. 
The distance function (9.3) can also being defined for almost-connections. 

The general properties of almost connections are being explored in another separated work. 
Some of these properties are based on straightforward generalizations of the quantities associated 
with Koszul connections. For instance, the generalization of the curvature tensor is: 

R(X, Y, Z) := L D t L DyZ - L D f L D t Z - L D [jt) Y]Z, X,Y,Z e TD. (11.2) 

The notion of auto-parallel curve of an almost-connection is not so easy to handle for almost- 
connections, since in general there will be initial conditions such that the auto-parallel curve 
goes out from the sub-bundle D. Indeed, for auto-parallel curves whose initial velocity vector 
is not on TD, the trajectory goes out from D after any finite time. 
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